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Abstract 

We investigate the effect of strong electron-electron repulsion on the electron- 
phonon interaction from a Fermi-liquid point of view. In particular we show 
that the strong interaction is responsible for vertex corrections, which are 
strongly dependent on the Vpq/uj ratio, where vp is the Fermi velocity and 
q and u are the transferred momentum and frequency respectively. These 
corrections generically lead to a strong suppression of the effective coupling 
between quasiparticles mediated by a single phonon exchange in the VFq/oj S> 
1 limit. However, such effect is not present when VFq/uj -C 1. Analyzing 
the Landau stability criterion, which involves the effective interactions in the 
dynamical limit, we show that a sizable electron-phonon interaction can push 
the system towards a phase separation instability. 

A detailed analysis is then carried out using a slave-boson approach for 
the infinite-U three-band Hubbard model describing the basic structure of a 
CuC>2 plane in copper oxides. In the presence of a coupling between the lo- 
cal hole density and a dispersionless optical phonon, we explicitly confirm the 
strong dependence of the hole-phonon coupling on the transferred momentum 
versus frequency ratio. We also find that the exchange of phonons leads to 
an unstable phase with negative compressibility already at small values of the 
bare hole-phonon coupling. Close to the unstable region, we detect Cooper 
instabilities both in s- and d-wave channels supporting a possible connec- 
tion between phase separation and superconductivity in strongly correlated 
systems. 

PACS:71.27.+a, 63.20.Kr, 74.72.-h, 74.25.Kc 



I. INTRODUCTION 

The problem of the interplay between the electron-phonon (e-ph) coupling and the 
(strong) electron-electron (e-e) interaction is an interesting problem, which still lacks a 
complete understanding. At present this topic is very hot since various facts indicate that 
the lattice can play a non negligible role in both the superconducting copper oxides and the 
fullerenes. As far as these latter materials are concerned, the strength of the interaction 
is still a matter of a debate, but there is a rather large agreement on the prominent role 
of the lattice. On the other hand, as far as copper oxides are concerned, it is generally 
recognized that the e-e interaction is very large in these systems, but it is the relevance of 
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the lattice that is questioned. However various groups 0-0] claim that there are both 
experimental evidences and theoretical arguments supporting a prominent role of the lattice 
in the cuprates. Moreover recent optical experiments in the mid-infrared frequency region 
indicate the presence of polaronic effects J| for the very lightly doped compounds, which 
are known to be strongly correlated systems. 

Various issues can be addressed in investigating the role of the lattice in the presence of a 
strong e-e interaction. In particular two questions can be raised, which are relevant, both on 
a general ground as well as in the framework of high temperature superconductivity. The first 
question concerns the possibility of having a large phonon-mediated effective e-e coupling 
due to the strong mass enhancement (m*/m ^> 1) occurring in a strongly correlated Fermi 
liquid: after all the dimensionless coupling A = j 2 uq (7 is the usual bare electron-phonon 
coupling and uq is the free-electron density of states) could grow very large because of 
the density of states renormalization arising from the mass enhancement A — > A = 7V*, 
with v* = (m*/m)z/ . This would have crucial consequences both on transport properties 
as well as on the Cooper pair formation. Moreover, a large A would favor the formation of 
polarons. 

The second question regards the possible occurrence of instabilities in the electronic gas. 
Previous studies of single- and multi-band Hubbard models in the strong-coupling (U ^> t) 
limit have revealed a strong tendency of these systems to undergo Phase Separation (PS) and 
Charge-Density- Wave (CDW) instabilities as soon as short range interactions are introduced. 
This occurs irrespectively to the magnetic (e.g. nearest neighbor Heisenberg coupling) ||- 
P2|l or the Coulombic (e.g. nearest neighbor repulsion) |I3|1 - [fL5] nature of the short-range 



interaction. In the framework of interest here one can ask whether also an e-ph coupling 
can destabilize the electron gas. The analysis of the stability with respect to PS of a given 
model is particularly relevant since the investigation carried out in models displaying PS 
also showed that superconductivity takes place close to the instability region ||, §|, |12|| - 
]rj| (the Cooper pairing occurring as a precursor of PS due to the attraction eventually 
driving the system to the PS instability). 

We anticipate here the answers to the two above questions, which will be analyzed in 
detail in the rest of the paper. 

As far as the first problem is concerned, a general analysis performed within a standard 
Fermi-liquid scheme with phonons coupled with the electron density reveals a dependence of 
the effective phonon-mediated e-e interaction on the ratio between the transferred momen- 
tum q and the frequency u, because of the vertex corrections generated by the e-e interaction. 
This dependence can be particularly strong in the presence of a large quasiparticle mass en- 
hancement. As a result, when vpqjuj ^> 1, we find that the screening of the quasiparticles 
strongly enhances the vertex corrections. These corrections tend to suppress the effective 
interaction so that the resulting static coupling A is small. This is in agreement with previ- 
ous calculations |T7|| performed in a specific model, where, however, the e-ph coupling arises 
from the ion-position dependence of the hopping matrix elements (the so-called covalent 
e-ph coupling). On the contrary, in the opposite limit vpq/uj <C 1 the intraband quasipar- 
ticle screening is ineffective, the physics being dominated by high energy (e.g. interband 
charge-transfer (CT)) processes, and a large effective e-ph interaction results. 

The strong dependence of the vertex corrections on the v^q/cu ratio renders the analysis 
of the effects of the e-ph coupling particularly delicate, since different physical quantities may 
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involve different vfq/oj regimes. In particular the e-ph scattering in transport properties is 
dominated by low-energy-high-momentum processes, which suffer a strong suppression due 
to the large vertex corrections. As pointed out in Ref. P7|, this will reduce the relevance 



of the e-ph scattering particularly at low doping, where the e-e interaction effects are more 
relevant. Similar conclusions have also been recently drawn |jTB[ from the analysis of a 
single-band Hubbard model with electrons coupled with an optical dispersionless phonon. 

On the other hand an affirmative answer can be given to the second question concerning 
the occurrence of instabilities in the electron gas. In fact, while the behaviour of A described 
above is generally true far from an instability, it will be shown that a sufficiently large e-ph 
coupling can instead produce a PS instability |L9|]. This reflects the fact that the Landau 
stability criterion for the symmetric Landau parameter, Fq > —1, required for a positive 
compressibility, involves the total (e-e and e-ph mediated) interaction in the dynamical limit 
Vpq/uj — > 0. At the point where F S = —1 the compressibility diverges and nearby, even 
the total static scattering amplitude (within a RPA resummation of the phonon-mediated 
interaction) is large and negative. Also this analysis can be carried out on a rather general 
ground. However, the quantitative determination of the instability conditions, being related 
to the subtle interplay between various interactions, depends on the couplings involved 
and must rely on the treatment of a specific model. For this purpose we shall consider a 
three-band Hubbard model describing the holes in the Cu02 planes of the high T c copper 
oxides. A coupling between the local hole density and a dispersionless optical phonon will be 
specifically considered. To deal with the electronic correlations in the strong-coupling limit 
we use the standard slave-boson technique within a 1/N expansion. The analysis of this 
model explicitly detects the presence of an instability region where the compressibility of the 
fermion gas diverges and then becomes negative. Near the region where the system becomes 
unstable, superconducting instabilities are found by averaging over the Fermi surface the 
interaction amplitudes in the Cooper channel. 

The plan of the paper is as follows. In Section II we address the two above questions 
working in the general framework of the Landau Fermi-liquid theory. The three-band Hub- 
bard model with the electrons coupled with a dispersionless optical phonon is introduced in 
Section III. The instabilities are investigated in Section IV. Final remarks are contained in 
Section V. 



II. PHONON-MEDIATED EFFECTIVE INTERACTION: A FERMI-LIQUID 

DISCUSSION 

One of the basic concepts of Landau Fermi-liquid theory is the idea of mass renormaliza- 
tion. This idea is also crucially present in the most common treatments of the interacting 
Fermi systems. In particular it is naturally introduced in the Gutzwiller treatment of the 
Hubbard model [ETI] and is at the basis of the Mott-Hubbard transition in this model [22 



This very same basic concept recurs in other popular techniques like, e.g. the slave-boson 
technique J29|- |2j. 



In the presence of a strong interaction the mass renormalization can be very large, 
m*/m ^> 1, and results in an enhancement of the quasiparticle density of states v* = 
(m*/m)u . Then, as already discussed in the previous Section, the natural question on 
the consequences of this renormalization on the effective phonon mediated e-e interaction is 
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whether the bare (free electron) effective e-e coupling A = 7 2 z/o grows into A = 7 2 z/* as an 
effect of the mass enhancement. 

Without a significant loss of generality we address the above question discussing the 
case of an optical phonon coupled to the local density of electrons by a constant coupling 
7 |24[. In the presence of a (possibly large) e-e interaction one has to worry about the 



e-ph vertex corrections involving the e-e interaction, for which no Migdal theorem can be 
applied. The problem is conveniently cast in the language of the standard Fermi-liquid theory 
|23| , by using the two relations connecting the density vertex A(q,u) and the wavefunction 



renormalization z w in the dynamic and static limits 

z w A(q = 0,ou^O) = 1 (1) 
z w A(q^0,uj = 0) =— ^ - (2) 

Fq = 2u*Y U) is the standard Landau parameter and Y u is the dynamic (q = 0, to — > 0) 
effective e-e scattering amplitude between the quasiparticles. 

We first consider the effective e-e interaction arising from a single phonon exchange 
at lowest order in r y 2 . Then the vertex corrections can not include phonon processes and 
the dynamic Landau scattering amplitude in Eq.fl2|) is due to the e-e interaction only. To 
explicitly keep memory of this limitation we append a suffix "e" to Y u and to any quantity 
not involving phononic processes. Thus F s(e) = 2z/T* . The relations (@) and (g) are exact 
Ward identities, which must be satisfied irrespective of the details of the e-e interaction and 
show a drastic difference between the dynamic (q = 0,u — > 0) and static (q — > 0,u = 0) 
limits. Whenever the exchange of a phonon takes place, the vertex corrections must be 
included leading to a different behaviour of the effective interaction in the two limits. The 
effective dimensionless e-e interaction mediated by a single phonon exchange reads 

SI$(q,u>) = z/7 2 4 2 A e2 (^) (3) 

where w(q) is the phonon dispersion. Here the presence of z e w indicates that we are consid- 
ering the effective interaction between quasiparticles and A expresses the difference of the 
phonon coupling to the quasiparticles with respect to particles |23| . 



The effects of the strong oo — q dependence of the electronic density vertex A, in Y e ff(q } u) 
[see Eq. (|3|)) can be made apparent in the small q and to limits, where the relations (|l|) and 
(H) can be used. Then one obtains 

i/TSSfe-O^-O) =-tV, — -0 (4) 



UJ 

V* 



where n e = -, — — -r is the compressibility of the Fermi liquid in the absence of coupling 

with the lattice. It should be noted that, in the present case of an optical phonon, the 
phonon propagator tends to -1 in both limits, Eqs.(^j) and (|5[). 

The difference between the dynamic and static case can be dramatic in the case of a 
Fermi liquid with a large mass enhancement m*/m ^> 1, but with a negligible compressibility 
renormalization (n e ~ v$) [pHj. In fact, in the case under consideration, Fq is proportional 



to the quasiparticle density of states, u* = (m*/m) u ^> i/ , and one has Fq ^> 1. Then 
Eqs.(0) and <E) read 



i/ItS(g = 0,wO) =-7 2 ( — U (6) 

*-pp/i 

yq — r u,uj — v) — — J — 

V 

(7) 

so that the effective one-phonon-mediated e-e interaction is large (~ m*/m) in the dynamic 
limit and small (~ m/m*) in the static one. 

The strong uj — q dependence in Eqs.(||) and (0) concerns the small-g and small-a; limits. 
This result relies on quite general arguments, whereas the case of finite g's and tu's will 
need the analysis of a specific model and will be discussed in the context of the three-band 
Hubbard model in the next sections. Our expectation, which will be confirmed by the 
analysis in Section IV, is that the product z w A will be roughly of order one (dynamical 
limit) all over the region outside the particle-hole continuum, while it will strongly deviate 
from unity in the region of the particle-hole continuum, where screening processes take place. 
As a consequence the e-ph coupling (and the e-e interaction mediated by phonons) will be 
depressed by the e-e interaction in all processes involving small energy and large momenta 
(as in the low energy lifetime and transport). 

The above conclusion on the irrelevance of the e-ph coupling is based on a lower order 
analysis in r y 2 . It evidently contrasts with the fact that the Landau stability criterion 
Fq > — 1 involves the full e-e interaction in the dynamical limit, where the e-ph coupling is 
not depressed by the "pure" e-e vertex corrections. Within the same limits discussed above 
(lowest order in 7 2 ) we have 

Fo = (K ~ 7 2 ) (8) 

indicating that a sizable 7 2 can indeed lead to F S < — 1. (Note that m*/m ^> 1 requires 
a large bare repulsion in units of the bare Fermi energy, as in the single- or multi-band 
Hubbard model near the metal-insulator transition, but this does not imply a large T^). The 
lowest order analysis showing the depression of the e-ph coupling in the low energy processes 
maintains its full validity with respect to the inclusion of higher order terms provided Fq in 
Eq.([8]) is still of order m*/m. On the contrary, near the instability condition F S = —1, the 
phonon contributions to the vertex cannot be neglected and the e-ph interaction is relevant 
even in the static limit. To clarify this point we extend our analysis considering a bubble 
resummation of the one-phonon processes, in the framework of a standard RPA approach 
for the quasiparticles. 
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As a first step we introduce an effective dynamic e-e scattering amplitude T^(q, uj) with no 
phonon processes included in it. Here and below, the subscript uj indicates that no intraband 
screening processes are present in the considered quantity. The q — > and uj — > limit of 
r^(g, uj) reproduces the above Landau dynamic scattering amplitude of the Fermi liquid 
in the absence of e-ph coupling. 

The second ingredient of our analysis is the intraband Lindhardt polarization bubble 



/ (E k+ i) - f [E k _L, 
k E k+ i - bj k _± - uj 

where / (E k ) is the Fermi function and E k is the quasiparticle band. Here we have included 
the spin multiplicity in the definition of H(q, uj). This bubble can be dressed with all possible 
intraband and interband purely electronic screening processes within a RPA resummation 
leading to 

W ' } l + TZ(q,co)U(q,u) 1 ; 

Finally we consider the self-energy corrections to the phonon propagator arising from the 
purely electronic processes. Once both the interband and the intraband screening of the 
quasiparticles is taken into account one obtains a phonon propagator with a self-energy 
correction given by g 2 Il 



uj 2 — UJ 2 



(q) (l- 7 2 n e (g,a/ 
uj 2 (q) 



, ,2 , ,2( n \ l+r u ,(g,aQII(g,aO 



(11) 



l+rg,(g,w)n(9,w) 

Notice that the quantity 

r u (q, uj) = rt(q, uj) + Tf(q = 0,uj^0) = T^q, uj) - 7 2 (12) 

has been introduced. Taking the limit q —>■ first and then uj — > 0, V becomes the Landau 
dynamic scattering amplitude in the presence of an e-ph coupling. 

We now can evaluate the total scattering amplitude splitting the contribution of the 
purely electronic processes and the contribution of the processes also involving phonons 

^ T(au)= SIZ(q,u) 1 ^ 2 (g) 7 2 fl3l 

l + TZ{q,u)U{q,u) [1 + r*(q, uj)U(q, uj)] 2 ^ _ u 2, s ^(^)n(^) 

It is important to notice that vertex corrections have been included in the second term 
by specializing Eqs.(^) and (||) into z e w K e = [1 + T^(q, u)H(q, uj)}^ 1 ■ Nevertheless, now 
the attraction mediated by this second term can be large even in the static limit if the 
denominator of the phonon propagator becomes small. In fact, in the uj — > limit one has 
for the attractive part of Eq. flT5| ) 
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(14) 



[1 + TZ(q, u)U(q, u)] 1 + r w (g, 0)1% 0) ' 

From this expression one can see that the condition for a diverging static scattering ampli- 
tude, eventually leading to an instability of the Fermi liquid, results in the condition on the 
Landau dynamic scattering amplitude 

1 + f LJ (q -> 0,lu = 0)1% -> 0,u = 0) = (15) 

equivalent to the usual 1 + Fq = 0, Eq.(|8|), since in this limit H(q,u) reduces to 2zA 

This result shows that a Fermi liquid can indeed be destabilized by the coupling of the 
quasiparticles to the lattice even in the presence of a very strong e-e repulsion. It must be 
noticed that this instability does not imply the full phonon softening, since the determination 
of the renormalized phonon frequency &(q) requires the value of 

1 + r\,(g,u;)ri(g,^) 

l + TZ(q,u)U(q,u) 1 ' 

in a range of frequencies, where the numerator does not vanish. The instability appears 
instead as an overdamping of the zero sound driven by the e-ph mediated attraction. 

The quantitative determination of the needed strength of the e-ph coupling in order to 
have an instability must rely on the analysis of a specific model. In fact, as it is apparent 
in its definition, Eq. fll2D , the strength and the sign of depend on the balance between 
T^(q,u) and 7 2 . Moreover, as we shall explicitly see in the model treated in the following 
sections, r^(g, u) in turn results from a cancellation between the strong bare repulsion and 
the strong interband screening. How much is left from this cancellation depends on the 
specific model one is dealing with. 

III. THE MODEL 

The model we consider is represented by the following Hamiltonian 

ia ia <ij>a 

+ J2 (tppijPtxPjay + h - C +U d J2 n di^dil + ^oJ2 A l A i 
<ij>a i i 

- (A + 4) [Gd (n* - (n di )) + G p {n pi - (n pi ))] . (17) 



where e d and e® are the Cu and O energy levels respectively, = ±t p d is the Cu-0 hybridiza- 
tion, t pp ij = ±t pp is the nearest neighbor 0-0 hybridization (for the choice of the orbital 
phases and the related choice of the sign of the hopping constants, see, e.g., Ref. |2B|| ). Ud is 
the on site repulsion between holes on copper sites. Starting from a Cu(3c/ 10 )-O(2p 6 ) vacuum 
state, holes on copper d x 2_ y 2 orbitals or on oxygen p x or p y orbitals at site i are created by 
the d\ and p\ operators respectively, = J2a d\ a di a is the total density per cell of holes on 
copper, while n pi = J2a,a=x, y PiaaPi,<r,a is the total density per cell of holes on oxygen. The 
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boson creation operator A\ creates a dispersionless phonon with frequency uoq coupled to the 
local density of copper and oxygen holes by the coupling constants Gd and G p respectively 
pTfl . These couplings can arise from the dependence of the interionic Coulombic repulsion 
on the relative position as well as on the hole occupation. An estimate can be obtained 
when only the nearest neighbor Coulomb interaction V(\r{ — rj\)ndin P j is considered. Then, 
a first order expansion on the ion displacement and on the charge fluctuations leads to an 
e-ph coupling of the form appearing in the Hamiltonian fllTD, with 



G d oc 



dV 

dri 



( n pj) i 



G p oc 



R 



dV 



(ridi) 



(18) 



R 



where R is the equilibrium distance between nn ions. Therefore, two different constants, 
Gd and G p , have been introduced owing to the different average occupations of the d and 
p orbitals and to the difference in ionic masses entering the standard normalization factors 
implicitly included in the definition of the G's. However, we do not want to stress very much 
in the present context the relevance of the specific form of the e-ph coupling appearing in 
Eq.fllBD p8| . In fact we do not address the specific aspects (like, e.g., symmetries and 
strengths) of the phonons in the cuprates, our issue being the understanding of the general 
properties of an Holstein [IE|| phonon in a strongly correlated system. 
A similar three-band Hubbard model was considered in Ref. 



171, where, however, an 



intersite "covalent" e-ph coupling was considered arising from the ion-position dependence 
of the hopping integrals. 

Since our investigation concerns the interplay between strong interactions and the 
phonons we take the Ud — > oo limit. In a standard way we handle the no-double-occupancy 
constraint on copper sites by means of the slave-boson technique p9|- fl32| . Therefore, 
after performing the usual substitution d\ a — > d i(T bi, d^ — > b\di a the constraint becomes 
J2a dicrdia + b\bi = 1. To equip the model with a formal small expansion parameter, we 
introduce a standard large-N expansion ||30|| , where the spin index a runs from 1 to N. The 
constraint is relaxed to assume the form J^a d\ a di a + b\bi = — and the suitable rescaling of 
the hopping t p d — > t p d/VN must, in this model, be joined by the similar rescaling of the 
hole-phonon coupling G — > G/yN in order to compensate for the presence of N fermionic 
degrees of freedom. Once these transformations are carried out the partition function of the 
final model can be written functional integral 



Z = J D P l a Dp aa DdlDd a Dtf DbDXDADA^- So Sdr , 



ST n't 9dia I \' J dpiaa jdbi ^dAi 
(j aa=x,y " ' L " 



+ E 



i\i b\bi - 



N 



(19) 
+ H, (20) 



// Y, d l* d ™ ( '/ + /A <) +4 53 i >>.*.<> - ■■■> 



i,a 
t 



t,cr,ot=x,y 



rrj [fPi,o-,a: Pi,a,—x Pi,a,y Pi,cr,—yj di a bi + C.( 

V-'" 1,(7 

tpp [pl+x,cr,a; {Pi,cr,-y Pi,cr,y + Pi+2x,cr,y Pi+2x,cr,—y) + C.C.J 



(A + 4) [G d (n di - (n dl )) + G 



(21) 
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where a local Lagrange multiplier field A has been introduced to implement the local con- 
straint forbidding the double occupancy on copper. 

At the mean-field (N = oo) level, the model of Eqs. (|T9|)-(|2~T|) is equivalent to the standard, 
purely electronic three-band Hubbard model without coupling to the phonons, which has 
been widely considered in the literature [Q, pEfl . In fact, at mean field level no role is played 
by the phonons because our electron-lattice coupling depends on the difference between 
the local and the average density and this difference naturally vanishes in the mean-field 



approximation [|34| . The model displays a T=0 Fermi-liquid behaviour for any finite doping 
5, where 5 is the deviation from half-filling, when one hole per cell is present in the system. 
In the Fermi-liquid case the mean-field value of the slave-boson field 6 multiplicatively 
renormalizes the hopping thus enhancing the effective mass of the quasiparticles (bo/y/N < 
Moreover, at this level the single particle self-energy does not introduce a finite 



quasiparticle lifetime. Then, in this model the single particle Green function of the physical 
fermions at iV = 2 has a quasiparticle pole with a finite residue given by the square of the 
mean-field value of the slave-boson field b^. 

On the other hand, at half-filling the system becomes an insulator if the bare charge- 
transfer energy difference e p — e° d — At pp is larger than a critical value ranging from 3.35t pd 



when t pp = [33] to a smaller value ~ 2.5t„d when t pp = 0.5t p d ||26||. In the insulating 



phase b vanishes leading to an infinite quasiparticle mass (m*/m = oo) and to a vanishing 
quasiparticle spectral weight |35] . 



In order to get new physical effects from the presence of the coupling with the phonons, 
one needs to consider the fluctuations of the bosonic fields. Since only a particular combina- 
tion a = (A^+A) I (2\/~N) of the phonon fields A and A^ is coupled to the fermions, it is conve- 
nient to use the field a and to integrate out the orthogonal combination a = (A— A 1 ) / (2\/ r N). 
Then the quadratic action for the boson field a reads 



H pbm = N y £ f £±^4a i , (22) 

n,i U 

where we have transformed the imaginary time into Matsubara frequencies. Moreover, 



working in the radial gauge [pi] , the phase of the field bi = v iVr* exp(— i(j>) is gauged 
away and only the modulus field is kept, while A, acquires a time dependence A; — > 
Aj + d T (pi. Thus one can define a three-component field A^ 1 = (<5r, 5X, a) where the time- 
and space-dependent components are the fluctuating part of the boson fields r { = r (1 + Sri), 
Xi = —iXo + 5Xi and aj. Writing the Hamiltonian of coupled fermions and bosons as H = 
HMF+Hb os +H- mt , where Hmf is the mean-field Hamiltonian quadratic in the fermionic fields, 
if bos is the purely bosonic part, also including the terms with the a, r and A bosons appearing 
in the action (^) and in ifphom Eq.(^). H int contains the fermion-boson interaction terms. 
More explicitly, Fourier transforming to the momentum space, the bosonic part reads 

H hos = Nj2^(q)B^(q)A»(-q) 

without explicitly indicating the frequency dependence for the sake of simplicity. The form 
of Eqs . (|20| ) - (|22| ) allows to determine the matrix B^ ,v ', whose elements are all zero except for 
B 1 ' 1 = r 2 A , B 1 ' 2 = B 2 ' 1 = *r 2 , 5 3 - 3 = (u 2 n + oo 2 ) fu . 
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To simplify the notation in Hmf and in H mt , we introduce a three-component fermionic 
field ^kaa = (d ka .,ip xka ,ip yka ) describing in momentum space the fermions in the orbital 
basis. Then Hmf can compactly be written as Hmf = Y^kaap H^ F {k)^>\ aa ty ka p, where 



Hmf (k) 



-2r t pd sm(k x /2) 
\ -2r t pd sm(k y /2) 



-2r t pd sm(k x /2) -2r t pd sin(k y /2) \ 



J 



with f3 k = 2sin(k x /2) sin(k y /2) and e d = £ d + Ao being the mean-field-renormalized en- 
ergy level of the copper d orbitals. The above matrix can be put in diagonal form by 
a unitary transformation to the quasiparticle basis ^ k aa = J2i3U a /3(k)^ka/3 leading to 
Hmf = Ukaa E a (k)^/l aa ^/ ka . a , where E 12i3 (k) are the mean-field-renormalized quasipar- 
ticle bands (we choose the band 1 as the lowest one, so that the Fermi level lays in this band 
for filling 5 < 1). This formalism allows to write the boson- fermion interaction as 



k,q,a 



(23) 



k,g,a 



The fermion-component index has been dropped and the (3 x 3) boson-fermion interaction 
vertices A^ 1 in the orbital operator basis can be obtained from Eq. (|2l|) 



A 1 



-2r n t 



Ql'pd 




sin 



sin {hz^ \ 










0\ 




' -2G d 


A 2 = 





, A 3 = 


-2G p cosf 




{o o oy 




v 









-2G P cos f 



(24) 



while the quasiparticle vertices A^g (k, q) are defined as 

A M (k, q) = u(k+'^\ A^ (k, q) f/ f (k 

Introducing the fermionic bubbles coupled to the various bosons 

f (E a (k + $)) - f (E p (k - * 



(25) 



E 

k,a,/3 



E a (k+I) 



Ep{k 



lUJr 



A^(k,q)A^ a (k, -q). 



(26) 



one can define the boson propagator 



D^(q,u m ) = < A^(q,u m )A u (-q,-uj m ) > 
= N~ l {2B + Yi{q^ m ))-l 



(27) 



The factor 2 multiplying the boson matrix B arises from the fact that the bosonic fields in 
the presently used radial gauge are real. 
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The set of formal tools is then completed by the introduction of the effective scattering 
amplitude between the quasiparticles in the lowest band 



r(fc, k'; q,iu) = -J2 Mi -q) D"" (?, w) (k, q) . (28) 

Then the scattering amplitude in the Cooper channel is given by 

Y c (k } k'~ uj) = -J2 K (^-, -k + fc') D^(q = k-k\ iv)A» n k - */) (29) 

It should be noted that the boson propagators are of order 1/N while the occurrence of a 
bare fermionic bubble leads to a spin summation and is therefore associated with a factor N. 
Thus, in this 1/N approach, the quasiparticle scattering amplitudes are residual interactions 
of order 1/N. The matrix form of the static density-density correlation function at the leading 
order is 

P a p(q,u = 0) = — ^ < n aa (q)n P(jl (-q) > 

era' 

= P° p (q, u = 0) + N X°^(q, u = 0) D^(q, u = 0) X %(q, w = 0) (30) 

where 

P ap(q^) = -^J2< n aa{q)n^{-q) > (31) 
era' 

is the orbital bare density- density correlation function, and 



X° a „(q, ^) = ^ E < n ™{q) E *L 7 A£(A;, q^k+qv's >o (32) 



1 

r 

crcr' fc;7,<5 

where a = d,p x ,p y , and /i = 1, 2 and 3. 



IV. THE DYNAMICAL ANALYSIS 

We now analyze the interplay between lattice and holes in the strongly correlated system 



represented by the model of Eqs.(|T9|)-(|2~T|) focusing on the possible occurrence of instabilities. 

Two mechanisms for driving a system to an instability can in principle be devised. A 
first one requires the complete softening of a collective mode leading to a ground state with 
different symmetries. For instance, in our specific model, this mechanism could imply the 
softening of the phonon leading to a structural transition. To investigate this mechanism a 
dynamical analysis of the collective modes of the system will be carried out in Subsection 
A. 

However, as we generally found out in Section II, an instability can also occur when the 
Landau criterion for stability Fq > —1 is violated (and usually a PS will then take place). 
In this case, the instability shows up as an overdamping of the zero-sound mode, the other 
collective modes remaining massive, and is a result of the delicate unbalance between the 
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repulsive and the attractive forces present in the system. Subsections B and C are devoted to 
this particular aspect, focusing on the role of the various screening processes in determining 
the instability. 

The result of the analyses performed in Subsections A, B and C is that the model of 
Eqs.([T^)-(|2T|) undergoes an instability driven by the second of the above mechanisms, and, 



in general, no phonon softening is detected. A similar phenomenon |T3|- [T5|] occurs in the 
three-band Hubbard model in the presence of a nearest neighbor Coulombic repulsion. In 
this latter case it is the attraction due to copper-oxygen CT fluctuations, which determines 
a PS instability without being accompanied by any collective-mode softening. 



A. The collective modes 

The collective modes appear as poles of the density susceptibilities. However an impor- 
tant remark is that the poles of the density susceptibilities \ coincide with the poles of the 
boson propagator D^ u (q,u), which therefore contain all of the relevant informations on the 
collective modes. In particular, to find the dispersion of these modes to = ou(q), one has to 
solve the equation 

det (2B + n) = 0. (33) 

Quite generally the charge-density fluctuations in the three-band model can be decomposed 
into the fluctuations of a non-conserved field, n d — n p , coupled to a fluctuating conserved 
field, rid + n P - The mode describing the propagating total density fluctuations is the zero- 
sound mode. This mode is massless, since in the model of Eqs.(|2l1|)-(|22|), no long-range 
interaction among the holes is included, which would otherwise transform this mode into a 
massive plasmon. In the presence of an overall attraction (0 > F$ > —1), the zero-sound 
mode enters into the particle-hole continuum and becomes damped. On the other hand the 
p-d CT fluctuations, being described by a non-conserved field, will contribute to the CT 
"optical-like" mode. Of course these fluctuations are dynamically coupled to the phonon 
mode. 

In order to simplify our analysis we first consider the q = limit because in this limit all 
intraband fermionic bubbles vanish thereby leaving the boson propagator in a much more 
treatable form. Physically this is due to the fact that in the q = limit the total density 
field p(q) = n d (q) + n p (q) = (l/y/W s ) Ek ^ +S(7 A + (fc, q)^ k -%a with 



2 



/l 

cosf | (34) 
\ cos f 



decouples from the dynamics because of particle conservation. To allow for a simpler analytic 
treatment we neglect in this section the direct oxygen-oxygen overlap (t pp = 0). Thus the 
boson-quasiparticle vertices simplify to the form 



A? 2 (M = 0) = ^f^ 



( A 

i 

\2G 
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M 3 (k,q = 0) 



Rk 



( 


\2G 



(35) 



where A = — Ed is the mean-field-renormalized atomic level difference and G = G p — Gd- 
Then the following expression for the polarization bubble results 



W v (q = Q,w) 



-?r 2 
z ' 



/ A 2 zA 2AG\ 

iA -1 2iG I(uj) 
\2AG 2iG AG 2 J 



where 



E 

k 



N x ^t 1 Rk(Rl 



oj 



Using Eq.(^Tj) one obtains the following expression for D 1 



D- L (q = 0,u) 



2Nr 2 Q 



t A - A 2 / i(l - A/) 

i(l - A/) 7 
v -2AJG -21TG 



-2AIG 
-2UG 
-AIG 2 + ^ 



and the determinant can be evaluated as 



det^ZT 1 (0,u;) 



2r 2 



uj q — OJ 

LJ 



;i-/(2A-A ))-4G 2 / 



(36) 



The analytical treatment can then proceed further introducing the small r 2 , approximation. 
This is justified at low doping deep inside the region of parameters where the system is 
insulating at half-filling. We recall that the model we are considering at leading order in 
1/N has a vanishing r 2 , approaching the insulating regime. Then it is convenient to rewrite 
I(oj) as 



with 



at(u)) 



A 



A 2 



UJ" 



1 



E 



a(u)r 2 



v iv s k 

so that the determinant assumes the form 



AniV, Rk (Ri 



UJ" 



(37) 



UJ Q — u 



(A - A ) 



+ A r 2 -AG 2 u + K - ^) ( 2A " A o) a(u) 



(38) 



By expanding a at zeroth order in r 2 



a (a;) 



16tL 74 



1 



7 2 A 2 - cj 2 



(39) 
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with 72 = (l/i\g J2kll = 1/2 + 2/tt 2 ~ 0.702 and 7 4 = (l/N a ) Eklt = 5/8 + 4/vr 2 ~ 1.03 
one obtains the expression of the determinant at first order in r$. It is immediate to recognize 
that, taking the r — > limit, one obtains two collective modes 



(A-X f 



0J„ 



(40) 
(41) 



corresponding to the excitonic charge-transfer mode (n p — n& fluctuations) and to the bare 
phonon mode. At this leading level in the 1/N expansion the two modes do not mix. On 
the other hand the equation det D~ x = can be solved giving 



uj 



± 
1 

~ 2 
+ rl 



+ v exc + r C (u) A 

2 



2 2 



+ \ r 2 [2C(u) ( 



2 2 
^exc — ^0 



16G 2 lu 



(u 2 + 



U exc) =•= ~~ ^excl 



C(u7) A ± 



Ao 



^0 _ ^ezcl 



C (E7) 



a; 2 + 8G 2 a; 



(42) 



where the short notation C (to) = (2A — Ao) oi (u = uJ) was introduced and the result was 
expanded at first order in r 2 ,. ZU assumes the values u exc or c<Jo corresponding to the zeroth 
order term in the expansion of to in r 2 . 

Two cases can be distinguished in Eq . (|42|) : oj exc > u , the more physical one, and 
uj exc < ujq which can only occur close to the metal-charge-transfer-insulator transition, where 

m 



the exciton mode completely softens [15 



In the former case one has 



UJ 



UJ 



u 2 exc + rfXo [C (uj exc ) + 



AG 2 uj q 



r 2 A 



AG 2 u 



whereas in the latter case the modes are 



W 2 = ^Ixc + r o% C {u e 



AG 2 u 



UJ 



4G 2 ^ 



0'MJ 2 2 
^0 — 10 exc 



(43) 
(44) 

(45) 
(46) 



It should be noted that in both "level repulsion" occurs so that the higher energy 

mode is pushed at even higher energy, while the lower energy one is made softer by the 
reciprocal interaction. However, it is important to stress the fact that this q = analysis 
does not reveal any complete softening of the collective modes in the limit of vanishing r 2 , 
where the energies of the modes reduce to their "bare" values uj exc and ujq. This indicates 
that any occurrence of q = instabilities in the small r region of the system is not due to 
a softening of the modes. Of course the above analysis was confined to the q = case and, 
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therefore, it does not allow to draw any conclusion on the occurrence of a collective mode 
softening at a finite wavelength or at large r . The investigation of this possibility requires a 
numerical analysis which was indeed performed. The result is that the static susceptibilities 
diverge first at q = and the softening of the modes at finite q only takes place inside an 
unstable region, characterized by F£ < — 1. This instability occurs as a consequence of the 
repulsion vs attraction unbalance generally presented in Section II, which we now elucidate 
for the specific three-band model under consideration. 



B. The Landau Fermi-liquid analysis 



In order to investigate the instability conditions of the model we now turn to a Fermi- 
liquid analysis transposing our 1/N calculation to the standard formalism presented in Sec- 
tion II and identifying the Landau amplitudes at leading order in 1/N. In particular, starting 
from the definition (^) of the singlet effective scattering amplitude one can introduce the 
standard Landau amplitudes 

r u (k, k') = - lim lim £ An (A/, -q) (q, u) A" n (k, q) . 

To carry out an analytical treatment we again assume t pp = 0. Then, since the model only 
includes a direct copper-oxygen transfer integral, T^k, k') depends on k and k' only via jp. 
and 7^'. Therefore, by taking the quasiparticles at the Fermi surface (where 7fc=constant= 
jp), only the "zeroth" harmonic would be non zero and given by 



r w = - hm lim £ (kp, -q) D»» (q, u) {k F , q) . 



(47) 



To explicitly get the expression of we first calculate 



2Nr 2 



1-: 



2R kp 



2R k 



V "2 \pdui F + G p vi F 



J 



( 



x 



a — i — ia(A — A ) 

— i — ia(A — Ao) A — (A — A ) 2 a 
2r ah(uj) 



2r a/i(u;) 
-2zrg(A - X )ah(uj) 



-2ir 2 (A - X )ah(u) rgu; (1 + 4r 2 a/i(u;)) / (uq - 
( K F -* 2 \ 

2R kp 

x ? ; ^ +A 

2R k p 



V "2 [Gdvlp + G p v 2 kF 



(48) 



/ 



with a = a(u) = [1 - (2A - A + 8r 2 V) I(u)} and h(u) = w G 2 /(^o ~ This S ives 

r(w) = r + a(a;)ri(u;), (49) 

with 
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ATr 



-L ( 1 + A_ ) {2Rkp _ 2A + Ao) ^ 

$rg \ R kp J uj 2 



UJ 



2 y G dU 2 kF + G>£ F 



(50) 



1 






2? 







(i?A 



2A + A ) 



+ (G d - G p ) (G d u 2 kF + G p v 2 k 

UJq — UJ N 



(51) 



A simple expression for can be obtained taking the u — > limit of the above formulas 
in the small r limit. It is worth noting that, whereas a(uj) is always finite, both T and r x 
grow very large when r — > 0. This is reminiscent of the infinite U we started with. However 
a straightforward expansion of Eqs.(^) and (|5l|) shows that the leading terms of order 1/rg 
(which characterize the model even in the absence of the e-ph interaction) cancel and only 
finite contributions are left. Therefore in this model a finite effective scattering amplitude is 
the result of a cancellation between a very large bare repulsion and a very large attraction 
due to interband screening, to which a finite phonon contribution is added. 
The resulting dynamical scattering amplitude is, for —>■ 



AT, 



4A 



(Ao - A) L 



2t 



A 2 



pd 



K^, 



(XqG p — AGd) 
2Xoujo 



(52) 



The first term represents the purely electronic (interband screened) repulsive contribution 
to [|HJ , whereas the second term is the attraction, which arises from the e-ph interaction. 
In particular it would coincide with the second term in the r.h.s. of Eq . (|12D , if the simple 
case G p = Gd = g is considered and the identification 7 2 = g 2 /uj is made PBfl . One can easily 
check that the e-ph vertex corrections in Eq.(|52D (i.e. the factors A a Ap with a, (3 — p,d, 
which multiply G a Gp/ujQ) can be written as 



d£° 



G„p(kF, oj 



0) 



(53) 



where G qp (k,uj) 



[UJ 



Et) 1 is the quasiparticle Green function and the derivatives are 



taken at fixed doping 5. One can then discover that in a multiband model the vertex 
corrections are still different from zero in the dynamical limit (cf. Eq.(^) with z w — 1 
because we are considering quasiparticles) and in fact can act to enhance the vertices. For 
instance, at 5 = + , we obtain A£ = Ao/(A — A), which reduces to unity only in the 
A —>■ oo limit, i.e. e° p — e° d —>■ oo. To obtain the vertex correction in the static limit we need 

to evaluate ^-o-G'Mkp, uj = 0) at fixed chemical potential fl3SH . We obtain 



A? 



i + Nv*r e „ 



(54) 



where is the first term in the r.h.s. of Eq. (p2\j . Eq.(Q) shows the dramatic difference 
between the dynamic and static limits. As discussed in Section II, this is a generic feature of 
strongly correlated systems, which one has to take into account since it can strongly affect 



the relevance of the e-ph coupling [40] 
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C. The phase diagram 



The previous analysis of the restricted model with t pp = provides valuable informations 
on the possible occurrence of instabilities of the Fermi liquid. In particular, having found 
r^, one can determine the Landau parameter F$ = Nh>*r w . Since the criterion for a finite 
positive compressibility is Fq > —1, and since v* oc \ 3> 1 in the small doping region, 



'0 



close to the insulating region (e^ — e d > 3.35t pd ) the stability criterion reduces to > 0. 



Taking into account the self-consistency condition Ao = 4^72/ A one is lead to the following 
inequality 

(X G P - AG d ) 2 < U 2 pd X u { ^ - 1 + ^ ] (55) 

Even for a weak e-ph coupling, this condition can easily be violated (see below). 

A clearer understanding of the stability condition can be obtained by considering the 
phonon coupled separately with the copper or the oxygen local density. In the case of 
Gd = the condition (^) becomes 



whereas for G p = one has 



^0 A \72 




Then one can easily recognize that, in the positive doping case with Ao > A, the coupling 
of the phonon with the oxygen holes G p is more effective in driving the instability. The role 
of A and A is interchanged at negative doping so that A < A and in this latter case the 
coupling Gd plays a major role in rendering the Fermi liquid unstable. Notice, however, that 
the model itself is not particle-hole symmetric so that the expressions (56) and ( |57D do not 
interchange in going from negative to positive doping. 

It is worth noting that, no matter how small the e-ph coupling is or how large the phonon 
frequency is, the instability will always take place by increasing the large bare charge-transfer 
difference (e p — e d )/t p d- This is due to the fact that, for large (e p — e d )/t p( i one has A ^> A 
for positive doping or Ao A for negative doping. Then it is possible to make the r.h.s. of 
Eqs.(|56|) or (|57|) smaller than the l.h.s.. 

The estimates (jj6]),([5"7|) were derived for the three-band Hubbard model with t pp = 
close to the insulating regime, where r ~ 0. To investigate the finite-doping regime and 
the model with finite oxygen-oxygen hopping we explicitly carried out a leading order 1 /N 
analysis of the static g-dependent density-density response function x(q, u; = 0) =< (n p (q) + 
n d(<l))( n p(—<l)+ n d(—<l)) >, which gives the compressibility of the system once the q — > limit 
is taken. This response function is obtained as a linear combination of the orbital density- 
density correlation functions in Eq . (|30|) . A divergent xi*!^ — 0) signals an instability of 
the Fermi liquid occurring at a wavevector q. 

We performed the calculations for many values of the parameters, and specifically for 
various phonon frequencies u and various e-ph couplings G p and Gd, and we particularly 
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explored the range of parameters, which could be relevant for the copper oxides. According to 
recent estimates [Q , [Q we considered various values of u Q ranging in the interval between 
0.01 and 0.08eV, particularly focusing on the 0.02eV region, where the phonon density of 
states is highest. As far as the e-ph couplings are concerned, these quantities are not directly 
accessible in experiments and their theoretical calculation is not easy due to the strongly 
interacting nature of the high temperature superconductors and to the nesting occurring in 
their Fermi surface [ |45[| . This is why we explored many possible cases obtaining qualitatively 
similar results. A reasonable choice is to use both G p and Gd different from zero, because a 
coupling of the lattice with both oxygen and copper hole density is naturally expected. The 
values of Gd and G p were such that a reasonable effective coupling A < 1 results far enough 
from the instability region. Moreover the choice G p > Gd appears reasonable in the case 
of intermediate hole doping according to the observations reported in the paragraph after 
Eq.(|l7D: i) G p is proportional to (n d ), while Gd is proportional to (n p ), with {rid) usually 
being larger than (n p ); ii) a factor is included in the definition of G p , whereas an 

analogous factor 1 / a/ M Cu enters in Gd further justifying the assumption G p > Gd, since 
Mcu > Mo- Both theoretical calculations and experiments in superconducting cuprates 
|J support the above qualitative arguments. The fact that oxygen holes are more strongly 
coupled to the lattice than holes on copper has relevant physical consequences in the light 
of the above observation that, for positive doping, G p is more effective than Gd in driving 
the instability. 

The phase diagrams in the {e p — e° d )lt p d vs 5 plane resulting from the analysis of x(q, w = 
0) are shown in Figs. 1 and 2 for typical parameter sets. The instability line indicates where 
x(q — > 0,uj = 0) diverges and it delimits a region of negative compressibility. The most 
impressive consequence of the coupling with the phonons is that a large unstable region 
appears at large doping or at large bare charge-transfer gap in both diagrams. It is worth 
noting again that the instability is not related to any collective mode softening. It can 
rather be interpreted as due to the attraction arising from the e-ph interaction eventually 
resulting in the overdamping of the zero-sound mode. Approaching the instability, the zero- 
sound mode first enters the continuum and gets Landau-damped. Eventually it becomes 
overdamped and the instability takes place. 

It is important to stress that the instability of the phase diagrams in Figs. 1 and 2 is an 
electronic instability even if it is driven by the coupling with the lattice. Moreover we find 
that the instability of x(q, uj = 0) first occurs at q — > and is therefore a signature of a long- 
wavelength, static, thermodynamical PS region to be identified by a Maxwell construction. 
To make this point more explicit, we report in Fig. 3 the Cooper-channel static scattering 
amplitude between quasiparticles on the Fermi surface T c (k F , k' F , u = 0) (see Eq.(^) and 
cf. also the to = limit of Eq.(|13"D) for a system with finite t pp (= 0.2t p d) and at various 
dopings. According to the observation already reported at the beginning of Subsection A, 
the scattering amplitude carries information on the instability because a divergency in this 
quantity can only be due to a divergency in the boson propagator also entering the expression 
for xiQ^ — 0), with q = k F — k' F . Thus a divergent scattering amplitude directly signals a 
divergent xiq^ = 0). As it can be seen in Fig.3, the divergency can also occur at sizable 
g's inside the unstable region. This can have relevant physical consequences, once a long- 
range Coulomb force is considered in the model. This latter interaction would stabilize the 
system in the regions of the phase diagram where x(<L uj = 0) shows divergencies at low-g. 
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This would prevent the occurrence of PS, but would leave open the possibility of finite-g 
instabilities: Most probably a sizable unstable region would survive to the introduction of a 
long-range Coulomb repulsion in the phase diagram parts where x(g, to = 0) was diverging 
at sizable g's, i.e. the divergency of the density-density response function would first occur 
at a finite q leading to the formation of incommensurate CDW. A similar phenomenon 



was first suggested for a different model in Refs. [fL3|l - [ 15| and seems to be confirmed by a 



static analysis of a pseudospin model in Ref. Of course the above argument misses the 
dynamical aspects of the problem, and should be taken as purely indicative of a possible 
scenario. One should also take into account that superconductivity can take place near the 
instability line where Fq = —1 and act to stabilize the system (i.e. superconductivity could 
partially or fully preempt the instability region). 

Few comments on the difference between the two phase diagrams of Figs.l and 2 are 
in order. It is apparent that the prominent effect of the direct oxygen-oxygen hopping is a 
stabilization at low doping and large e p — e d . In order to clarify this point we performed 
a calculation of the compressibility both in the presence and in the absence of t pp in the 
limit of very large (e p — e d ) / t p d and small positive doping. Due to the slave-boson- mean- 
field band renormalization the lowest quasiparticle band closely approaches the bottom of 
the intermediate band, mostly of oxygen character. This latter band is a flat non-bonding 
oxygen band in the absence of t pp and acquires a dispersion 4t pp in the presence of direct 
oxygen-oxygen hopping if the mixing with the copper d x 2_ y 2 orbitals is neglected. In the 
(e p — e d ) /t p d — > oo and 5 — > + one has a vanishing distance between the renormalized d 
atomic level and the bottom of the pure-oxygen band A = e p — 4t pp — Ed — > 0. In this limit 
the calculation greatly simplifies and, for t pp ^ the compressibility per spin assumes the 
form 

(58) 

t U ■ (59) 

b pp \ L pp 

This result is suggestive of a decoupling between the p and the d levels occurring when 
(£p — £°d)/tpd — ► oo: the compressibility for the interacting system of (1 + + ) holes in the 
mixed p-d quasiparticle band beares resemblance with the compressibility of a system of 5 = 
+ non-interacting holes in the pure-oxygen band. In this latter system the compressibility 
coincides with the density of states, which is of order 1 jt pp and is independent on doping in 
two dimensions, whereas it vanishes with doping in the d = 3 case. This is precisely the result 
reported in Eqs.(^8|) and (|59"D respectively, once one remembers that, in the (£ p — e d )/t p d ~^ oo 
and 5 — > + limits, A is vanishingly small. 

Eqs.(|58D and ( |59f) also express the fact that in the large (e p — e d )/t pd and low-doping limit, 
the compressibility of the system diverges when t pp tends to zero (even for the three-band 
Hubbard model in the absence of e-ph coupling). Since K e = u*/(l + iWT^), this result 
can be rephrased as ~ t pp . Thus, when t pp — > and the quasiparticle density of states 
v* is large, even a very small attraction leads to a negative T u , immediately resulting in a 
violation of the Landau stability condition F£ = N^T^ < —1. The three-band Hubbard 
model with t pp = and e p — e d = oo case is strictly analogous to the U = oo single-band 
case [[jjj . This clarifies why, the region with infinite compressibility at low doping and large 
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Ep — e° d in the phase diagram of the system of Fig.l (t pp = 0), is more extended than the one 
for the system of Fig. 2. 



D. The Cooper instability 

Finally we investigated the possibility of Cooper pairing. As already pointed out, super- 
conductivity can appear as a precursor of the PS instability. In particular a simple argument 
can be put forward fl3fl , fl4| suggesting that a region of large compressibility is a good can- 



didate in order to find superconductivity. In the general language of Landau Fermi-liquid 
theory the compressibility can be written as k = ^0"/ (1 + Fq)- Then, if the compress- 
ibility gets large without being accompanied by a large mass enhancement, (as it happens 
near the instability line), this means that the denominator is becoming small indicating a 
negative F$ (indeed F$ = — 1 at the instability line). Although this parameter is related 
to the quasiparticle scattering in the particle-hole channel, it seems reasonable to expect 
attraction in the particle-particle channel as well. According with this plausibility argument 
and according to the previous experience in other strongly interacting models, we therefore 
investigated the Fermi surface average of the particle-particle scattering amplitude defined 
in Eq.(g) 

_/ / dkdk'5 (E(k) -fi)S (E(k') - fi) gi{k)T c (k, k'; u = 0) gi (k') 
l ~ JdkS(E(k) - fi)gi(k) 2 ' 

with g si (k) = cos(k x ) + cos(k y ), g S2 {k) = (cos(k x ) - cos(k y )) 2 , g dl {k) = cos(k x ) - cos(k y ) 
and gd 2 {k) = sin(k x ) sm(k y ) projecting the interaction onto the s-wave and d-wave channels. 
(Notice that A; > means attraction). Whereas the couplings A^ are found to be generally 
attractive near (and inside) the unstable region, we find s-wave Cooper instabilities only very 
close to the instability line. The results are tabulated at various doping concentrations for 
the case with t pp = and (e° — £°)/t P d = 3.3i p d in Table 1 and for the case with t pp = 0.2tpd 
and (s® — e®)/t pd = 4:.9t p d in Table 2. In both cases the phonon frequency and the e-ph 
couplings are u = 0.02t pd , and G p = 0.15^, and Gd = 0.1t pd , respectively. With the set 
of parameters related to Table 1 the critical doping for the occurrence of the instability is 
5 C = 0.21, whereas for the set of Table 2 the instability is at 5 C = 0.23. Our analysis indicates 
the sure existence of d-wave pairing in sizable regions of model (^) in the limit of large-N 
near the instability, whereas the occurrence of s-wave pairing takes place in a much narrower 
region [J41]J . However, it should be emphasized that the presence of a s-wave static Cooper 
instability only in a narrow region, by no means excludes the possibility of having s-wave su- 
perconductivity in a much larger area of our phase diagram. Only an appropriate Eliashberg 
dynamical analysis can allow to draw a conclusion, specially in the light of the considera- 
tions on the strong frequency dependence of the vertex corrections discussed in Section II. 
Of course the same applies to the attraction in the d-wave channels, which could also be 
greatly favored by dynamical effects. The full Eliashberg dynamical analysis unfortunately 
involves a difficult complete treatment of the momentum and frequency dependence of the 
effective scattering amplitude, which is beyond the scope of the present paper. Nevertheless 
valuable informations can be gained from an analysis of the e-ph vertex. The aim of this 
analysis is to provide useful quantitative indications on the behaviour of the e-ph interac- 
tion as a function of both momentum and imaginary frequency in various parts of the phase 
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diagram as a preliminary step towards the implementation of the full Eliashberg treatment 
of the effective phonon-mediated e-e scattering amplitude. Also in this case we checked that 
our results do not qualitatively change by varying the model parameters. In particular we 
choose the values of the various quantities in order to allow for a possible connection with 
the experimentally known features of the copper oxides. We choose the value of the bare 
atomic level difference — e d = 4.9t p d so that a CT optical gap in the insulating phase 
of about 2t p d would result. Assuming the usual value of t p d =1.3eV for the copper-oxygen 
hopping integral, this latter value turns out to be 2.6eV, and is not much larger than the 
experimentally known values of the CT gap in the insulating phase of the superconducting 
copper oxides (e.g. in Ref. j|2| the reported value for La 2 Cu04 is 2.3eV). As values of the 
phonon frequency we take 0.02^ ~0.026 eV and 0.08t p d ~0.104 eV. The first value was 
chosen because in this frequency region the experimental phononic density of states is large, 
whereas the larger value was used to extend the region where dynamical effects can take 
place before the phonon frequency cuts off the effective interaction. The value of the e-ph 



couplings was chosen accordingly in order to obtain reasonable values of A [45 



We first report in Fig. 5 the interaction (slave-boson) renormalized vertex between 
phonons and quasiparticles without phonon processes included in it (see Fig. 4) as a function 
of the Matsubara frequency. Two of these quantities can be joined to a bare phonon prop- 
agator to give the one-phonon effective scattering amplitude specializing to the three-band 
Hubbard model the general expression (|3|). 

On the other hand, joining two of the above vertices with a fully renormalized phonon 
propagator leads to the full phonon-mediated scattering amplitude in the Cooper channel. 
This quantity was represented by the second term of Eq. (|i~3]) in Section II and is reported 
in Fig. 6 for the three-band Hubbard model as a function of the Matsubara frequency. The 
momenta of the external fermions correspond to the Cooper channel with k and k! on the 
Fermi surface. 

It is clear that the results reported in Figs. 5 and 6 quantitatively confirm the general 
qualitative analysis of Section II. Fig. 5 displays the strong dependence of the one-phonon 
vertex from the momentum vs frequency ratio. A typical bandwidth of about Q.2t pd was 
chosen, giving a Fermi surface density of states ~ 13/t pd . We only report the results cor- 
responding to a doping 5 = 0.225 since no significant dependence of this vertex from the 
doping is detected, so that the closeness to the instability is immaterial in this quantity. 
On the other hand, it is apparent that a rapid increase occurs in the phonon-quasiparticle 
vertex, when the frequency becomes larger than vpq. This is most evident at low momenta, 
where the phonon-quasiparticle vertex G e g increases fast by more than one order of mag- 
nitude, from the small (see inset in Fig. 5) value at zero frequency up to a sizable value at 
large-frequency. This latter is much larger because at frequency larger than the bandwidth 
for quasiparticle-quasihole excitations (vpq at small q in an isotropic system) the intraband 
electronic screening is ineffective (cf. Eq. ([[])) J46j. 

Analyzing in detail the contributions coming from the vertex corrections, we discovered 
that in this "high" -energy region the interband screening acts to enhance the vertex at small 
g's by a factor of about 1.5 with respect to the bare vertex (Go — 0.2). The bare value is 
only recovered at very high energy (much larger than the interband CT energy A), and/or 
large q. 

More dramatic are the momentum and doping dependences of the attractive scattering 



21 



amplitude in Fig. 6. In the stable region (5 = 0.2), a maximum in the attraction is observed 
at finite frequency and low momentum. This maximum is a result of the balance between the 
energy dependence of the vertex correction and of the boson propagator (see the inset). The 
interaction at low momenta is by far larger than the attraction at large momenta. The same 
behavior as a function of momenta holds near the unstable region of doping, (5 = 0.225), but 
at low energy and small momentum the attraction is not suppressed by vertex corrections, 
because in the vicinity to the instability line the attraction reaches its maximum at ou — 0. 

Finally Figs. 7 and 8 report the total effective scattering amplitude in the Cooper channel 
T c (k F , k' F , uj n ) (continuous lines) resulting from the sum of the attractive amplitude (dotted 
lines) of Fig. 6 and the repulsion arising from purely electronic scattering processes (dashed 
lines) as a function of Matsubara frequencies at small (a) and large q (b). 

Being the calculations performed close to the instability, we see that the e-ph interaction 
gives a large contribution to the total interaction at small q, leading to an overall attraction 
at low frequencies. On the contrary, at large q the full interaction is only governed by the 
pure e-e repulsion and no phonon energy scale is detectable in its frequency dependence. 
Fig. 8 reports the same quantities of Fig. 7, but for a larger phonon frequency (consequently 
G p has been increased, G p = 0.3t p d, to drive the system unstable at reasonable doping 
5 C = 0.185). Qualitatively the results of Fig.8 show the same behaviour depicted in Fig. 7, 
but it is evident that the phonon-mediated attraction extends over a larger frequency range. 
We also found that in this latter case the attraction persists up to larger momenta. At large 
q (Fig.8b) the results for T c \k F , k F ,u n ) are still very similar to those of Fig. 7b. 

Finally Fig.9 reports T c (k F , k' F ,u n ) as well as its attractive and purely electronic repul- 
sive components at 5 = 0.2. Due to the depressing effect of the vertex at low frequencies, the 
phonon contributes very little to T c (k F , k' F , u n ) even at small q. Nevertheless, as it is also 
apparent from Fig.6, at frequency of the order of u the attractive part contributes sizably 
to the total interaction. This clearly indicates that at small q, contrary to the static attrac- 
tion, which is strongly suppressed by the vertex corrections, the finite frequency attraction 
persists in regions of doping that are far from the instability line. 

The analysis summarized in Figs. 5-9 is of obvious pertinence in a complete Eliashberg 
treatment of the superconductivity problem. In particular it is evident that the huge en- 
hancement of the attractive part of the scattering amplitude near the instability line can be 
responsible for large critical temperatures despite the small e-ph coupling. The closeness to 
a PS instability appears therefore as a favorable condition in order to obtain high temper- 
ature superconductivity from a phonon-mediated attraction similarly to what suggested in 
the context of purely electronic pairing mechanisms. 

V. CONCLUSIONS 

In the present work we investigated the e-ph interaction in the presence of a strong local 
repulsion within a general Landau Fermi-liquid framework. Using standard Ward identities, 
we pointed out the strong dependence of the e-ph vertex on the momentum vs frequency 
ratio in a strongly interacting system displaying a large effective-mass enhancement, but not 
too large compressibility. In particular we showed that the dimensionless attractive quasi- 
particle scattering mediated by a single phonon exchange coupled to the electronic density 
is strongly suppressed by vertex corrections due to the e-e repulsion when v F q > oj. On the 
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contrary the dimensionless effective one-phonon-mediated attraction is strongly enhanced by 
the effective mass increase when vpq < to. These results stay valid within a random-phase 
resummation of phonon exchanges, provided the system is far from an instability This 
strong momentum vs frequency dependence of the e-ph coupling shows up in different ways 
in different physical quantities. In particular it was found in a three-band Hubbard model 
with intersite "covalent" e-ph coupling [|17]], that phonons contribute little to the quasipar- 
ticle scattering as far as transport properties are concerned. We expect a similar depression 



for the e-ph coupling considered in this paper [49]. This is so because transport properties 



involve low-frequency-high-momentum processes, where the large-momentum limit of the 
e-ph interaction is mostly involved. In this limit the e-ph interaction is strongly suppressed, 
particularly in the low-doping regime, where the effects of the strong e-e interaction play 
a major role. A different behaviour is expected in other quantities where large frequencies 
are more relevant. This is in agreement with recent calculations |18| performed in a single- 
band Hubbard model with an on-site "ionic" e-ph coupling, which show, indeed, that the 
Eliashberg function a 2 F{uj) determining superconductivity, is much less reduced than the 
analogous quantity determining the transport properties. 

These remarks are obviously relevant in any dynamical analysis of the pair formation in 
strongly interacting systems like, e.g., high temperature superconducting cuprates, fullerenes 
or Bai_ x K x Bi03. In particular, our work calls for a critical reanalysis of effective potential 
models in the Eliashberg approach of pairing || ]5(| , when a strong e-e interaction is present 
in the systems. Most of the work in this field has, in fact, been carried out using model 
potentials, where the momenta are averaged on the Fermi surface. This leads to a trivial 
momentum dependence (likely overestimating the role of large momenta), which misses the 
peculiar strong momentum vs frequency dependence revealed in the interacting systems. 

On the other hand our general analysis showed that the e-ph coupling can drive the 
strongly correlated system to an instability. This was specifically shown for the three-band 
Hubbard model in the U — > oo limit, where more specific features arise. In particular it was 
evidenced that the Fermi-liquid system has diverging long-wavelength density fluctuations 
in some regions of the — E®)/t pi i vs. 5 phase diagram, even with a reasonably small 
e-ph coupling. This is due to the phonon induced attraction leading to the violation of the 
Landau stability criterion Fq > — 1. In analogy with other cases Jl3|- [|l^] this appears as an 
overdamping of the zero-sound mode. Moreover it was found that superconducting pairing 
both in the s- and d-wave channels occurs close to the instability region. It is worth stressing 
that the detected Cooper instabilities occur in the presence of an infinite local repulsion U. 

These findings are of obvious theoretical relevance in the context of high temperature 
superconductivity. The general considerations of Section II show that one is not allowed to 
draw any naive conclusion on the extension of the superconducting region on the basis of a 
static analysis, which neglects the presence of the e-e interaction in the e-ph coupling, once 
strong correlations are present. However, it should be noticed that the results of Section II 
indicate that, if a system displays Cooper phenomenon in a static (u = 0) analysis including 
vertex correction, the pairing can only be favoured by the extension of the analysis to finite 
frequencies. We also like to point out that phonon corrections to the e-ph vertex, which 
are usually neglected according to the Migdal theorem [|7|, have recently been considered 
in Ref. [48]. This analysis showed that these corrections, which are not negligible if the 



bandwidth vs phonon frequency ratio is sizable, tend to suppress the vertex in the large 
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VFq/oj region, whereas they tend to enhance the e-ph coupling when vfq/uj is small. These 
corrections tend therefore to cooperate with the strong interaction effects discussed in the 
present paper. 

Our results suggest that, if high temperature superconductivity is driven by e-ph cou- 
pling, this is possibly due to the fact that the system is close to a PS instability, where the 
argument of the vertex suppression does not apply. If this is the case, superconductivity 
itself would possibly provide the system a stabilizing mechanism with respect to supercon- 
ductivity 0. 

By increasing the bare e-ph coupling, one could indeed get a sizable attractive \ s (u = 0) 
at small and intermediate doping if, for instance, the e-ph mediated attraction is evaluated 
at order 7 2 (a similar calculation has been carried out in Ref. |5l| for the single-band infinite- 
U Hubbard model with a covalent (intersite) e-ph coupling). However, at the same time, Fq 
(which is less affected by e-e vertex corrections), evaluated in the normal state, should be 
largely negative (F S <C —1) and one should worry about the stability of the system. The 
straightforward approach would be to allow for pairing and then evaluate |^ and then check 
for thermodynamic stability. This requires the evaluation of e-e screening effects in the 
superconducting state, where hopefully the (partial or full) opening of the gap will stabilize 
the system. In this paper we have kept the aptitude to evaluate the dimensionless couplings 
A; in the stable region where Fq > —1, leaving to a future work the analysis of the region 
with Fq < — 1 in the normal phase. 

The generality of the arguments presented in Section II and the qualitative independence 
of the results obtained in the framework of the three-band Hubbard model by varying the e- 
ph couplings G p and Gd and other parameters of the model (like, e.g., the phonon frequency 
ojq or the oxygen-oxygen hybridization t pp ), witness that our results are rather robust and 
should also depend little on the particular choice of the phonon considered in the model. The 
choice taken in the model (pj]) was dictated by simplicity, but we do not expect qualitative 
differences in the physics of more realistic models. Besides superconductivity, our results 
provide a possible explanation for other general features of the superconducting cuprates 
as well. In particular our results concerning the presence of incommensurate CDW would 
explain the formation of the superstructures, which are present in many copper oxides, if 
this phenomenon has an electronic origin. 

A relevant issue related to the work presented here concerns the possible formation of 
polarons in the three-band Hubbard model. The treatment of this long standing problem 
is beyond the scope of the paper, but it is worth emphasizing that the band narrowing 
occurring at low doping due to strong e-e repulsion could ease the polaron formation. 

Some insight on this can be gained from the analysis of Monte Carlo calculations |52]] 
showing that a polaronic regime may take place in a single-band model with one single 
electron coupled via a short-range interaction g to optical phonons. The polaronic regime 
occurs as soon as the e-ph coupling g exceeds a critical value g c . The numerical analysis 
shows that the critical value depends on both the electronic bandwidth (~ t in Ref. ||52|| ) and 
the phonon frequency (u ), and it can be represented by the condition g^/(tco ) = X c ~ 1 
in the adiabatic, t 3> uq, regime or by the condition g c /u>o = a c ~ 1 in the antiadiabatic, 
t <C uj Q , regime. 

The numerical Monte-Carlo analysis |52|] only concerns the case of one single electron in a 
lattice, whereas the presence of other interacting electrons may modify the above picture. As 
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discussed in Section II, the high sensitivity of the vertex corrections as a function of frequency 
and momentum renders the above conditions for polaron formation rather ambiguous, since 
it is not clear which frequency vs momentum regions dominate the e-e screening processes 
dressing the e-ph coupling g (cf. Section II). A quick inspection to the (cfcj(a, + aj)) 
correlation function shows that a rather natural guess is that the relevant frequencies are 
smaller than the phonon frequency u>o- Then, if the quasiparticle bandwidth t* is larger than 
co>o, the vertex corrections to the phonon propagator are largely to be considered in the large 
momentum (vpq/uj > 1 and Vpq ~ t*) and are large, thus leading to a suppression of the e-ph 
correlation. This would disfavor the establishment of a polaronic regime. On the other hand, 
as soon as the effective fermion bandwidth t* decreases (e.g. approaching the zero-doping 
insulating phase) the vertex corrections are in the high frequency region (t* m vpq < oj) 
and become small, so that the effective e-e coupling is given by A = g 2 /(ut*), which grows 
large. This increase likely leads to polaron formation. This considerations suggest that, if 
the e-e correlations do not lead to a fermion effective bandwidth smaller that the typical 
phonon frequencies, polaron are more difficult than in the non interacting system. On the 
contrary, if the closeness to an insulating phase reduces the effective bandwidth determining 
non-adiabatic conditions for the quasiparticles, then the polaronic regime is made easier: 
the e-e interaction itself reduces the electronic kinetic energy thus favoring the gain in 
lattice deformation energy related to the polaron formation. Notice that the lattice-driven 
instability of the Fermi liquid, which we showed to occur in the three-band Hubbard model 
for not too large e-ph coupling, does not appear to be related to the formation of polarons 
in the system. In fact the instability occurs at small g's (specifically at q = 0), while the 
formation of local polarons involves all gs. Moreover the polaron instability is associated to 
the rapid softening and hardening of the phonon mode, while our dynamical analysis has 
shown that the instability is a consequence of an overdamping of the zero-sound, the phonon 
remaining massive. 

A last remark can be made concerning the possible formation of bipolarons. In the limit of 
weak e-e repulsion (U ~ 0), and in the antiadiabatic regime, a large enough phonon-mediated 
attraction would lead to a violation of the Landau stability criterion Fq < —1 (since ~ 0, 
r^" < and v* is large) signaling the instability of the Fermi liquid. This criterion, instead, 
is not violated in the low-doping limit of the infinite-?/ three-band Hubbard model (at least 
in the more physical case with a finite oxygen-oxygen hybridization). This indicates that, 
even when the antiadiabatic conditions are realized (close to the insulating phase), the large 
repulsion is able to stabilize the system against bipolaron formation: if polarons are formed, 
the low frequency analysis indicates that they are stable against bipolaron formation. Of 
course this does not exclude the possibility af a dynamical binding of the polarons. 

A clear confirmation of our results concerning the presence of a lattice-driven instability 
comes from the numerical exact diagonalization of the three-band Hubbard model with e-ph 
coupling in a small cluster |53|] . This calculations show that a CDW occurs at finite doping 
(one hole in a cluster of four Cu atoms with surrounding oxygen) if the e-ph coupling does 
not exceeds a critical value. Above this value the system enters, instead, a small-polaron 
regime. It is quite appealing to associate the CDW in the small cluster to the unstable 
region in the infinite system considered in this paper. Of course, in the case investigated in 
the present paper by continuously increasing the doping, the instability first occurs at zero 
wavevectors, but this behaviour could hardly be detected in a finite small cluster, where 
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there are few available momenta and where the doping cannot be varied smoothly (the 
addition of a single hole already produces a doping around 0.25). 
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FIGURE CAPTIONS 

FIG. 1: Phase diagram (e p — e d )/t pd versus positive doping 5 for t pp = 0, G p = 0.15t p d, 
G d = O.ltpd, and u = 0.02t pd . 

FIG. 2: Phase diagram (e Q p —e° d ) /t pd versus positive doping 8 for t pp = 0.2t pd , G p = 0.15t pd , 
G d = 0.1t pd , and u = 0.02t pd . 

FIG. 3: Effective static quasiparticle scattering amplitude versus transferred momentum 
well outside (5 = 0.2, dotted line), slightly outside (S = 0.225, continuous line) and well 
inside (5 = 0.3, dashed line) the unstable region. The parameter are as in Fig. 2 with 
e p — e d = 4:.9t pd . The instability occurs at 5 C = 0.23. The scattered quasiparticles are on the 
Fermi surface and q is in the y direction, q = (0, q). 

FIG. 4: a) Diagrammatic representation of the renormalized quasiparticle-phonon ver- 
tex. The electronically screened vertex is reported in b), where the open dot is the bare 
phonon-quasiparticle vertex and the thick dashed line is the resummed slave-boson prop- 
agator diagrammatically represented in c). No phonon propagators are included in the 
resummation c). 

FIG. 5: Interaction-renormalized quasiparticle-phonon vertex at fixed low (q = 0.06, 
continuous line) and large momenta (q = 2.20, dashed line) as a function of the transferred 
Matsubara frequency oo. The parameters are as in Fig.3 with 5 = 0.225. The scattered 
quasiparticles are on the Fermi surface and q is in the y direction, q = (0,q). The inset 
shows an enlargement of the low frequency region. The quasiparticle density of states 
is v* ~ 13/tpd per cell. The vertical lines indicate the phonon frequency and half the 
quasiparticle bandwidth. 

FIG. 6: Effective dynamic full-phonon-mediated scattering amplitude at fixed low 
(q = 0.06, continuous line) and large momenta (q = 2.20, dashed line) versus transferred 
Matsubara frequency. The parameters are as in Fig.3. The calculations are performed well 
outside (5 = 0.2, diamonds) and slightly outside (5 = 0.225, crosses) the unstable region. 
The scattered quasiparticles are on the Fermi surface and q is in the y direction, q = (0, q). 
The quasiparticle density of states is v* ~ 13/tpd per cell. 

FIG. 7: (a) Effective dynamic quasiparticle total (continuous line), purely electronic 
repulsive (dashed line) and full-phonon (dotted line) scattering amplitudes in the Cooper 
channel at fixed low momentum (q = 0.06) as a function of the transferred Matsubara 
frequency. The parameters are as in Fig.3. The calculations are performed slightly outside 
(5 = 0.225) the unstable region. The scattered quasiparticles are on the Fermi surface and 
q is in the y direction, q = (0, q). The quasiparticle density of states is v* ~ 13/t pd per cell. 
The vertical lines indicate the phonon frequency and half the quasiparticle bandwidth, (b) 
Same as in (a) but at large momentum (q = 2.20). 

FIG. 8: (a) Effective dynamic quasiparticle total (continuous line), purely electronic 
repulsive (dashed line) and full-phonon (dotted line) scattering amplitudes in the Cooper 
channel at fixed low (q = 0.07) versus transferred Matsubara frequency. The calculations 
are performed at 5 = 0.18, close to the unstable region, occurring at S c = 0.185. The 
phonon frequency is uj = 0.08t pd , G p = 0.3t pd and G d = 0.1t pd . The scattered quasiparticles 
are on the Fermi surface and q = (0.053,0.03) (\q\ = 0.06). The quasiparticle density of 
states is v* ~ 9/t pd per cell. The vertical lines indicate the phonon frequency and half the 
quasiparticle bandwidth, (b) Same as in (a) but with large momenta q = (0,2.25). 

FIG. 9: Effective dynamic quasiparticle total (continuous line), purely electronic repul- 
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sive (dashed line) and full-phonon (dotted line) scattering amplitudes in the Cooper channel 
at fixed low momentum (q = 0.083) versus transferred Matsubara frequency. The phonon 
frequency is u = 0.02t pd , G p = 0.15t p d and Gd = 0.1t pd and the calculations are performed at 
5 = 0.2, far from the unstable region, occurring at 5 C = 0.225. The scattered quasiparticles 
are on the Fermi surface and q = (0.078,0.028) (\q\ = 0.083). The quasiparticle density of 
states is v* ~ 9/t p d per cell. The vertical lines indicate the phonon frequency and half the 
quasiparticle bandwidth. 
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TABLE CAPTIONS 

TABLE 1: Si- and c^-wave coupling constants for various doping and for the model with 
Bp — e d = 3.3t p d, t pp = 0, G p = O.lBtpd, Gd = 0.1t pd , and u = 0.02t pd . The instability line is 
at£ c = 0.21. 

TABLE 2: si-, s 2 -, d\- and o?2-wave coupling constants for various doping and for the 
model with e° — e d = A.9t p d, t pp = 0.2t pd , G p = 0.15t p d, Gd = 0.1t pd , and uj = 0.02t pd . The 
instability line is at 5 C = 0.23. 
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TABLES 



TABLE 1 



5 


0.15 


0.20 


0.208 


0.22 


A S1 


-0.75 


-0.5 


0.43 


-0.45 




-0.044 


0.2 


0.35 





TABLE 2 



5 


0.15 


0.20 


0.229 


Asi 


-0.58 


-0.48 


1.2 


A S2 


-0.55 


-0.65 


1.5 


Adj 


-0.063 


0.017 


2.1 


Ad 2 


-0.021 


0.012 


1.8 
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